(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 6.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     71328,       1939]
NotebookOptionsPosition[     65684,       1750]
NotebookOutlinePosition[     66027,       1765]
CellTagsIndexPosition[     65984,       1762]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[{
 RowBox[{
  RowBox[{"Off", "[", 
   RowBox[{"ClebschGordan", "::", "\"\<tri\>\""}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Off", "[", 
   RowBox[{"ClebschGordan", "::", "\"\<phy\>\""}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Off", "[", 
   RowBox[{"SixJSymbol", "::", "\"\<tri\>\""}], "]"}], ";"}]}], "Input",
 CellChangeTimes->{{3.456663504453125*^9, 3.45666350465625*^9}, {
  3.5407344053453703`*^9, 3.540734425171504*^9}}],

Cell[CellGroupData[{

Cell[TextData[{
 "\n",
 StyleBox["Two-photon line shape calculation for K-39 ", "Title"],
 StyleBox["D-2 line from 4S to 4D3/2", "Title",
  FontFamily->"Arial",
  FontVariations->{"Underline"->True}]
}], "Input",
 CellChangeTimes->{{3.551727958766342*^9, 3.5517280116923695`*^9}, {
   3.5517283803094535`*^9, 3.55172838042146*^9}, {3.5521680002246284`*^9, 
   3.552168000672654*^9}, 3.5527533840615234`*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"-", "39"}], " ", "D"}], "-", 
  RowBox[{
  "calculation", " ", "for", " ", "K", " ", "line", " ", "photon", " ", 
   "shape"}], "-", 
  RowBox[{"16", " ", "D3", " ", "from", " ", "line", " ", "S", " ", "to"}], 
  "+", "Two"}]], "Output",
 CellChangeTimes->{3.5521685414453125`*^9, 3.552753832883789*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Ouput Location ", "Subsubtitle",
 CellChangeTimes->{{3.5517158329247837`*^9, 3.551715854388012*^9}, 
   3.551716247196479*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{
   "oFile", "=", 
    "\"\<C:\\\\Users\\\\work\\\\Data\\\\output\\\\K39-D2-4D32.txt\>\""}], 
   ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{"output", " ", "file"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.5407383642928085`*^9, 3.5407383651238565`*^9}, 
   3.5407384082223215`*^9, {3.5407384849037075`*^9, 3.540738535987629*^9}, {
   3.54073857516387*^9, 3.5407386491461015`*^9}, {3.5407388655584793`*^9, 
   3.5407388798112946`*^9}, {3.5407453919391146`*^9, 3.5407453981011543`*^9}, 
   3.540745430736563*^9, {3.5409092806939077`*^9, 3.540909299115961*^9}, {
   3.5409139555021634`*^9, 3.5409139572030754`*^9}, {3.540916902879408*^9, 
   3.5409169748426695`*^9}, {3.5409174063863525`*^9, 
   3.5409174253559527`*^9}, {3.5413429497201304`*^9, 
   3.5413429616052127`*^9}, {3.5413431664128094`*^9, 3.541343179302828*^9}, {
   3.542556628837941*^9, 3.542556640226014*^9}, {3.54255670671364*^9, 
   3.5425567165729036`*^9}, {3.5425569015744896`*^9, 3.542556922993427*^9}, {
   3.5425654644126453`*^9, 3.542565470200283*^9}, {3.542570842449442*^9, 
   3.5425708492199287`*^9}, {3.542852782822707*^9, 3.5428527850223217`*^9}, {
   3.542854309978897*^9, 3.542854310868102*^9}, {3.54299342661618*^9, 
   3.5429934460288*^9}, {3.5429966016528063`*^9, 3.5429966440456142`*^9}, 
   3.5429966817575817`*^9, {3.5486039749375*^9, 3.5486040204375*^9}, {
   3.548604193265625*^9, 3.548604193390625*^9}, {3.549060869234375*^9, 
   3.54906087171875*^9}, {3.5492471925820312`*^9, 3.5492471933359375`*^9}, {
   3.5500073476933594`*^9, 3.5500073526621094`*^9}, {3.5512341310498047`*^9, 
   3.5512341453935547`*^9}, {3.5517158046171646`*^9, 3.551715826351408*^9}, 
   3.551716178698561*^9, {3.5517175755924587`*^9, 3.551717580352731*^9}, {
   3.5517284278251715`*^9, 3.551728427945178*^9}, {3.5521680129523563`*^9, 
   3.552168013395382*^9}, {3.552175399694336*^9, 3.552175401084961*^9}, {
   3.5527533777617188`*^9, 3.5527533778847656`*^9}}],

Cell[TextData[{
 StyleBox["Uncomment the code below to use ",
  FontFamily->"Times New Roman",
  FontWeight->"Bold",
  FontColor->RGBColor[1., 0., 0.]],
 StyleBox["experimental",
  FontFamily->"Times New Roman",
  FontWeight->"Bold",
  FontVariations->{"Underline"->True},
  FontColor->RGBColor[1., 0., 0.]],
 StyleBox[" reptition rate and offset frequency",
  FontFamily->"Times New Roman",
  FontWeight->"Bold",
  FontColor->RGBColor[1., 0., 0.]]
}], "Text",
 CellChangeTimes->{{3.5517161881951046`*^9, 3.551716263157392*^9}, {
  3.5517170355175686`*^9, 3.551717038884761*^9}, {3.5517277525975504`*^9, 
  3.5517277846923857`*^9}},
 FormatType->"TextForm"],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{
   "pathname", "=", "\"\<C:\\\\Users\\\\work\\\\Data\\\\2012-05-15\\\\\>\""}],
    ";", "\[IndentingNewLine]", 
   RowBox[{"filename", "=", "\"\<Rb-D2-line.dat\>\""}], ";", " ", 
   RowBox[{"(*", "datafile", "*)"}], "\[IndentingNewLine]", 
   RowBox[{"ExpSpecFile", "=", 
    RowBox[{"pathname", " ", "<>", " ", "filename"}]}], ";", 
   "\[IndentingNewLine]", 
   RowBox[{"dataTab", "=", " ", 
    RowBox[{"Import", "[", 
     RowBox[{"ExpSpecFile", ",", "\"\<Table\>\""}], "]"}]}], ";", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{"Form", " ", "repetition", " ", "rate", " ", 
     RowBox[{"list", ".", " ", "Extract"}], " ", "from", " ", "data", " ", 
     "file"}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{"fOffsetList", "=", 
    RowBox[{"{", "}"}]}], ";", "\[IndentingNewLine]", 
   RowBox[{"fRepList", "=", 
    RowBox[{"{", "}"}]}], ";", "\[IndentingNewLine]", 
   RowBox[{"n", "=", "1"}], ";", "\[IndentingNewLine]", 
   RowBox[{"While", "[", 
    RowBox[{
     RowBox[{"n", "\[LessEqual]", " ", 
      RowBox[{"Length", "[", "dataTab", "]"}]}], ",", "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"fRepList", "=", 
        RowBox[{"Append", "[", 
         RowBox[{"fRepList", ",", 
          RowBox[{
           RowBox[{"dataTab", "[", 
            RowBox[{"[", "n", "]"}], "]"}], "[", 
           RowBox[{"[", "2", "]"}], "]"}]}], "]"}]}], ";", " ", 
       RowBox[{"(*", 
        RowBox[{"\"\<2\>\"", " ", "=", " ", 
         RowBox[{
         "2", "nd", " ", "column", " ", "of", " ", "the", " ", "datafile", 
          " ", "contains", " ", "the", " ", "rep", " ", "rate"}]}], "*)"}], 
       "\[IndentingNewLine]", 
       RowBox[{"fOffsetList", "=", 
        RowBox[{"Append", "[", 
         RowBox[{"fOffsetList", ",", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "2"}], "*", 
            RowBox[{
             RowBox[{"dataTab", "[", 
              RowBox[{"[", "n", "]"}], "]"}], "[", 
             RowBox[{"[", "1", "]"}], "]"}]}], ")"}]}], "]"}]}], ";", " ", 
       RowBox[{"(*", " ", 
        RowBox[{
         RowBox[{"+", 
          RowBox[{"/", 
           RowBox[{"-", " ", "depends"}]}]}], " ", "on", " ", "locking", " ", 
         "electronics"}], "*)"}], "\[IndentingNewLine]", 
       RowBox[{"n", "++"}], ";"}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", 
   RowBox[{"Clear", "[", "n", "]"}], ";", "\[IndentingNewLine]", 
   RowBox[{"fRepMin", "=", 
    RowBox[{"Min", "[", "fRepList", "]"}]}], ";", "\[IndentingNewLine]", 
   RowBox[{"fRepMax", "=", 
    RowBox[{"Max", "[", "fRepList", "]"}]}], ";", "\[IndentingNewLine]", 
   RowBox[{"fRepAvg", "=", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"fRepMin", "+", "fRepMax"}], ")"}], "/", "2."}]}], ";"}], 
  "*)"}]], "Input",
 CellChangeTimes->{{3.5407383642928085`*^9, 3.5407383651238565`*^9}, 
   3.5407384082223215`*^9, {3.5407384849037075`*^9, 3.540738535987629*^9}, {
   3.54073857516387*^9, 3.5407386491461015`*^9}, {3.5407388655584793`*^9, 
   3.5407388798112946`*^9}, {3.5407453919391146`*^9, 3.5407453981011543`*^9}, 
   3.540745430736563*^9, {3.5409092806939077`*^9, 3.540909299115961*^9}, {
   3.5409139555021634`*^9, 3.5409139572030754`*^9}, {3.540916902879408*^9, 
   3.5409169748426695`*^9}, {3.5409174063863525`*^9, 
   3.5409174253559527`*^9}, {3.5413429497201304`*^9, 
   3.5413429616052127`*^9}, {3.5413431664128094`*^9, 3.541343179302828*^9}, {
   3.542556628837941*^9, 3.542556640226014*^9}, {3.54255670671364*^9, 
   3.5425567165729036`*^9}, {3.5425569015744896`*^9, 3.542556922993427*^9}, {
   3.5425654644126453`*^9, 3.542565470200283*^9}, {3.542570842449442*^9, 
   3.5425708492199287`*^9}, {3.542852782822707*^9, 3.5428527850223217`*^9}, {
   3.542854309978897*^9, 3.542854310868102*^9}, {3.54299342661618*^9, 
   3.5429934460288*^9}, {3.5429966016528063`*^9, 3.5429966440456142`*^9}, 
   3.5429966817575817`*^9, {3.5486039749375*^9, 3.5486040204375*^9}, {
   3.548604193265625*^9, 3.548604193390625*^9}, {3.549060869234375*^9, 
   3.54906087171875*^9}, {3.5492471925820312`*^9, 3.5492471933359375`*^9}, {
   3.5500073476933594`*^9, 3.5500073526621094`*^9}, {3.5512341310498047`*^9, 
   3.5512341453935547`*^9}, {3.5517158046171646`*^9, 3.551715826351408*^9}, 
   3.551716178698561*^9, {3.55171684287755*^9, 3.5517168751273947`*^9}, {
   3.5517169260393066`*^9, 3.5517169994325047`*^9}, {3.551727761577064*^9, 
   3.551727761648068*^9}}],

Cell[TextData[{
 StyleBox["Uncomment the code below to use ",
  FontFamily->"Times New Roman",
  FontWeight->"Bold",
  FontColor->RGBColor[1., 0., 0.]],
 StyleBox["theoretical",
  FontFamily->"Times New Roman",
  FontWeight->"Bold",
  FontVariations->{"Underline"->True},
  FontColor->RGBColor[1., 0., 0.]],
 StyleBox[" repitition rate and offset frequency",
  FontFamily->"Times New Roman",
  FontWeight->"Bold",
  FontColor->RGBColor[1., 0., 0.]],
 " "
}], "Text",
 CellChangeTimes->{{3.551716164954775*^9, 3.5517161669018865`*^9}, {
  3.55171627657716*^9, 3.551716297164337*^9}, {3.5517278230845823`*^9, 
  3.5517278563484845`*^9}},
 FormatType->"TextForm"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"f0", "=", 
    RowBox[{
     RowBox[{"-", "2"}], 
     RowBox[{"(", 
      RowBox[{"20.", " ", 
       SuperscriptBox["10", "6"]}], ")"}]}]}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"after", " ", 
     RowBox[{"SHG", ".", " ", "Add"}]}], " ", "+", 
    RowBox[{"/", "-"}]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"fRepMin", "=", "250000000."}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{"default", " ", "frep"}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"fRepRange", "=", "80."}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{"scanning", " ", "range"}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fRepMax", "=", 
   RowBox[{"fRepMin", "+", "fRepRange"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Delta]f", "=", "0.25"}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{"Hz", " ", "-", " ", 
    RowBox[{"step", " ", "size", " ", "in", " ", "rep", " ", "rate"}]}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fr", "=", "fRepMin"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fRepList", "=", 
   RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"fOffsetList", "=", 
    RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
   "Form", " ", "the", " ", "fRepList", " ", "and", " ", "fOffsetList"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"While", " ", "[", 
   RowBox[{
    RowBox[{"fr", "<", "fRepMax"}], ",", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"fRepList", "=", 
       RowBox[{"Append", "[", 
        RowBox[{"fRepList", ",", "fr"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"fOffsetList", "=", 
       RowBox[{"Append", "[", 
        RowBox[{"fOffsetList", ",", "f0"}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"fr", "+=", "\[Delta]f"}], ";"}], "}"}]}], "]"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"Average", " ", "fRep"}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fRepAvg", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"fRepMin", "+", "fRepMax"}], ")"}], "/", "2."}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.5517161466957307`*^9, 3.5517161587894225`*^9}, {
   3.5517163024016366`*^9, 3.5517163047197695`*^9}, {3.5517163366095934`*^9, 
   3.5517163419118967`*^9}, {3.55171653062969*^9, 3.5517167452519665`*^9}, {
   3.5517167950658154`*^9, 3.551716838278287*^9}, 3.5517170583998775`*^9, {
   3.5517286701030283`*^9, 3.551728670279039*^9}, 3.5521680173266068`*^9, {
   3.5521685349296875`*^9, 3.5521685350234375`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Physical Constants", "Subtitle",
 FontColor->RGBColor[0, 0.356863, 0]],

Cell["Speed of Light", "Subsubtitle"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"c", "=", "29979245800."}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"cm", "/", "sec"}], "*)"}]}]], "Input",
 CellChangeTimes->{3.551728122119686*^9}],

Cell[CellGroupData[{

Cell["Planck's Constant", "Subsubtitle"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"hBar", "=", 
    RowBox[{"1.054", " ", 
     SuperscriptBox["10", 
      RowBox[{"-", "27"}]]}]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"ergs", " ", "sec"}], "*)"}]}]], "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Boltzmann Constant", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"kB", "=", 
    RowBox[{"1.380658", " ", 
     SuperscriptBox["10", 
      RowBox[{"-", "16"}]]}]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"erg", "/", "K"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.4566596001875*^9, 3.456659625953125*^9}, {
  3.45665968075*^9, 3.456659686359375*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Atomic Parameters", "Subtitle",
 FontColor->RGBColor[0, 0.356863, 0]],

Cell["General Parameters", "Subtitle",
 CellChangeTimes->{{3.45666146840625*^9, 3.45666147296875*^9}}],

Cell["Nuclear Spin", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"i", "=", 
   RowBox[{"N", "[", 
    FractionBox["3", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.4730885615*^9, 3.473088561953125*^9}, 3.53632801824225*^9, {
   3.5425567297081876`*^9, 3.5425567322042036`*^9}, {3.5512806388623047`*^9, 
   3.5512806390654297`*^9}}],

Cell[CellGroupData[{

Cell["Transit Linewidth", "Subsubtitle",
 CellChangeTimes->{{3.45666389359375*^9, 3.456663897328125*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"\[Gamma]T", "=", 
   FractionBox["1", 
    RowBox[{"40.2", " ", 
     SuperscriptBox["10", 
      RowBox[{"-", "9"}]]}]]}], ";"}]], "Input",
 CellChangeTimes->{{3.456663726421875*^9, 3.4566637388125*^9}, {
   3.45666380203125*^9, 3.456663813921875*^9}, 3.4566922269375*^9, 
   3.5366835360235*^9, 3.538523997111776*^9, {3.5492471882460938`*^9, 
   3.549247189754883*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Polarization", "Subsubtitle",
 CellChangeTimes->{{3.456663909203125*^9, 3.456663910671875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"q1", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "1", ",", "0"}], "}"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.456663912765625*^9, 3.456663917453125*^9}, 
   3.456663957921875*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"q2", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "1", ",", "0"}], "}"}]}], ";"}]], "Input",
 CellChangeTimes->{3.4566639629375*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Ground State Parameters", "Subtitle",
 CellChangeTimes->{{3.456659852765625*^9, 3.456659860875*^9}}],

Cell["Ground State S", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
   3.45665970596875*^9, 3.456659711984375*^9}, {3.45665974584375*^9, 
   3.456659752203125*^9}, 3.45665980753125*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"SGnd", "=", 
   RowBox[{"N", "[", 
    FractionBox["1", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665975521875*^9, 3.456659757859375*^9}, 3.4566598095625*^9, {
   3.5425567378202395`*^9, 3.5425567406126575`*^9}}],

Cell[CellGroupData[{

Cell["Ground State J", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.45665974584375*^9, 
  3.456659752203125*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"jGnd", "=", 
   RowBox[{"N", "[", 
    FractionBox["1", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665975521875*^9, 3.456659757859375*^9}, {3.5425567479135046`*^9, 
  3.542556749863517*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Ground State F Values", "Subsubtitle",
 CellChangeTimes->{{3.45666290321875*^9, 3.456662908015625*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"fMinGnd", "=", 
   RowBox[{"Abs", "[", 
    RowBox[{"jGnd", "-", "i"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fMaxGnd", "=", 
   RowBox[{"jGnd", "+", "i"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.456662861609375*^9, 3.456662890515625*^9}, {
   3.54906118934375*^9, 3.549061201671875*^9}, {3.5490736629375*^9, 
   3.549073666609375*^9}, {3.5492478572304688`*^9, 3.549247859171875*^9}, {
   3.550007296041992*^9, 3.5500073055283203`*^9}, {3.5500073762246094`*^9, 
   3.5500073772871094`*^9}, 3.5512806623154297`*^9, {3.5512807609091797`*^9, 
   3.5512807639716797`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Ground State Hyperfine A", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
  3.45666013296875*^9, 3.456660152140625*^9}, {3.456660398109375*^9, 
  3.45666040421875*^9}, {3.456660756828125*^9, 3.45666075759375*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"aGnd", "=", 
    RowBox[{"SetPrecision", "[", 
     RowBox[{"230859860100.", ",", "20"}], "]"}]}], " ", ";"}], 
  RowBox[{"(*", " ", 
   RowBox[{"230.8598601", " ", "MHz"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
   3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
   3.456660078296875*^9, 3.456660092390625*^9}, {3.456660306609375*^9, 
   3.45666036596875*^9}, {3.456660407359375*^9, 3.45666043703125*^9}, {
   3.4566605581875*^9, 3.4566605719375*^9}, {3.456660680015625*^9, 
   3.456660705609375*^9}, {3.45666076115625*^9, 3.4566607615625*^9}, {
   3.47308857578125*^9, 3.473088584546875*^9}, {3.47309640846875*^9, 
   3.473096424140625*^9}, {3.536335670945375*^9, 3.53633567542975*^9}, {
   3.5512809188779297`*^9, 3.5512809360341797`*^9}, {3.5512810990966797`*^9, 
   3.5512811179873047`*^9}, 3.5512811900185547`*^9, {3.551715895359355*^9, 
   3.5517159071990323`*^9}, {3.551717675681184*^9, 3.5517176816565256`*^9}, 
   3.551728139904703*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Ground State Hyperfine B", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
   3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
   3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
   3.45666013296875*^9, 3.456660152140625*^9}, {3.456660398109375*^9, 
   3.45666040421875*^9}, 3.4566605775*^9, {3.45666076390625*^9, 
   3.45666076546875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"bGnd", "=", "0."}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
   3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
   3.456660078296875*^9, 3.456660092390625*^9}, {3.456660306609375*^9, 
   3.45666036596875*^9}, {3.456660407359375*^9, 3.45666043703125*^9}, {
   3.4566605581875*^9, 3.456660602203125*^9}, {3.45666063840625*^9, 
   3.45666063971875*^9}, {3.45666076803125*^9, 3.456660773828125*^9}, 
   3.5512811922998047`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 "Intermediate State Parameters ",
 Cell[BoxData[
  FormBox[
   RowBox[{"4", 
    SubscriptBox["P", 
     RowBox[{"3", "/", "2"}]]}], TraditionalForm]]]
}], "Subtitle",
 CellChangeTimes->{{3.456659852765625*^9, 3.456659887765625*^9}, {
  3.5512809799404297`*^9, 3.5512809843154297`*^9}, {3.5527533931523438`*^9, 
  3.5527533954726562`*^9}}],

Cell["Intermediate State S", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.45665974584375*^9, 
  3.456659752203125*^9}, {3.45665980753125*^9, 3.45665981840625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"SInt", "=", 
   RowBox[{"N", "[", 
    FractionBox["1", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665975521875*^9, 3.456659757859375*^9}, {3.4566598095625*^9, 
  3.45665982178125*^9}, {3.542556758911575*^9, 3.5425567611423893`*^9}}],

Cell[CellGroupData[{

Cell["Intermediate State J", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"jInt", "=", 
   RowBox[{"N", "[", 
    FractionBox["3", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665976928125*^9, 3.45665977578125*^9}, {3.473088588953125*^9, 
   3.47308858946875*^9}, 3.536333582476625*^9, {3.5425567651828156`*^9, 
   3.5425567679908333`*^9}, {3.5512809908779297`*^9, 
   3.5512809910029297`*^9}, {3.5527533995458984`*^9, 3.5527533996796875`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State F Values", "Subsubtitle",
 CellChangeTimes->{{3.45666290321875*^9, 3.456662923203125*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"fMinInt", "=", 
   RowBox[{"Abs", "[", 
    RowBox[{"jInt", "-", "i"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fMaxInt", "=", 
   RowBox[{"jInt", "+", "i"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.456662861609375*^9, 3.456662890515625*^9}, {
  3.456662927125*^9, 3.456662936640625*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State Linewidth", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Gamma]Int", "=", 
    RowBox[{"N", "[", 
     FractionBox["1", 
      RowBox[{"26.37", " ", 
       SuperscriptBox["10", 
        RowBox[{"-", "9"}]]}]], "]"}]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"Seconds", " ", "From", " ", "Budker"}], ",", "Kimball", ",", " ",
     "DeMille"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
  3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
  3.456660078296875*^9, 3.456660092390625*^9}, {3.53668310449225*^9, 
  3.536683106007875*^9}, {3.5512810147529297`*^9, 3.5512810175185547`*^9}, {
  3.5512811805185547`*^9, 3.5512811853779297`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State Center of Gravity Energy", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
  3.45666013296875*^9, 3.456660152140625*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"\[Nu]Int", "=", 
   RowBox[{"SetPrecision", "[", 
    RowBox[{
     RowBox[{"391.01617083", " ", 
      SuperscriptBox["10", "12"]}], ",", "20"}], "]"}]}], " ", 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Hz", " ", "from", " ", "Steck"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
   3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
   3.456660078296875*^9, 3.456660092390625*^9}, {3.456660306609375*^9, 
   3.45666036596875*^9}, {3.473088596515625*^9, 3.473088616875*^9}, {
   3.47308883840625*^9, 3.473088843484375*^9}, {3.473096433640625*^9, 
   3.473096439859375*^9}, {3.53633478836725*^9, 3.53633479311725*^9}, 
   3.536607319211*^9, {3.536607447226625*^9, 3.536607498007875*^9}, {
   3.53666646042975*^9, 3.5366664651485*^9}, {3.537893167558856*^9, 
   3.5378931893520937`*^9}, 3.5378943549667425`*^9, 3.5378974743522215`*^9, {
   3.5378985672711487`*^9, 3.53789856825395*^9}, {3.5378998431661897`*^9, 
   3.537899869374236*^9}, 3.5379017571978636`*^9, 3.538606320562795*^9, {
   3.5512810393310547`*^9, 3.5512810634560547`*^9}, 3.5512811714248047`*^9, {
   3.552753582833008*^9, 3.552753602392578*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State Hyperfine A", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
  3.45666013296875*^9, 3.456660152140625*^9}, {3.456660398109375*^9, 
  3.45666040421875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"aInt", "=", 
   RowBox[{"SetPrecision", "[", 
    RowBox[{
     RowBox[{"6.933", " ", 
      SuperscriptBox["10", "6"]}], ",", "20"}], "]"}]}], " ", ";", 
  RowBox[{"(*", 
   RowBox[{"Hz", " ", "From", " ", "Steck"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
  3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
  3.456660078296875*^9, 3.456660092390625*^9}, {3.456660306609375*^9, 
  3.45666036596875*^9}, {3.456660407359375*^9, 3.45666043703125*^9}, {
  3.4566605581875*^9, 3.4566605719375*^9}, {3.473088627078125*^9, 
  3.4730886284375*^9}, {3.473096443796875*^9, 3.473096448328125*^9}, {
  3.5363357313985*^9, 3.536335734726625*^9}, {3.536607079039125*^9, 
  3.53660709955475*^9}, {3.536666472757875*^9, 3.536666476257875*^9}, {
  3.536683146882875*^9, 3.5366831523985*^9}, {3.5372022874831247`*^9, 
  3.5372022919049997`*^9}, {3.5512811388466797`*^9, 3.5512811690185547`*^9}, {
  3.5517160806559534`*^9, 3.551716087167326*^9}, {3.5527536052421875`*^9, 
  3.5527536164472656`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State Hyperfine B", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
   3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
   3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
   3.45666013296875*^9, 3.456660152140625*^9}, {3.456660398109375*^9, 
   3.45666040421875*^9}, 3.4566605775*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"bInt", "=", 
   RowBox[{"SetPrecision", "[", 
    RowBox[{
     RowBox[{"2.786", " ", 
      SuperscriptBox["10", "6"]}], ",", "20"}], "]"}]}], " ", ";", 
  RowBox[{"(*", 
   RowBox[{"Hz", " ", "From", " ", "Steck"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
   3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
   3.456660078296875*^9, 3.456660092390625*^9}, {3.456660306609375*^9, 
   3.45666036596875*^9}, {3.456660407359375*^9, 3.45666043703125*^9}, {
   3.4566605581875*^9, 3.456660602203125*^9}, {3.45666063840625*^9, 
   3.45666063971875*^9}, 3.473088631296875*^9, {3.473096452546875*^9, 
   3.473096456171875*^9}, {3.536335742211*^9, 3.536335745961*^9}, {
   3.536607112711*^9, 3.536607116664125*^9}, 3.53666648311725*^9, {
   3.5366831642735*^9, 3.53668317030475*^9}, {3.5372023036862497`*^9, 
   3.5372023085924997`*^9}, {3.5512811663310547`*^9, 
   3.5512811667529297`*^9}, {3.55275362209375*^9, 3.5527536285117188`*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 "Final State Parameters ",
 Cell[BoxData[
  FormBox[
   RowBox[{"4", 
    SubscriptBox["D", 
     RowBox[{"3", "/", "2"}]]}], TraditionalForm]]]
}], "Subtitle",
 CellChangeTimes->{{3.456659852765625*^9, 3.456659887765625*^9}, {
  3.456660101109375*^9, 3.4566601018125*^9}, {3.5512819261904297`*^9, 
  3.5512819303154297`*^9}, {3.551717330572445*^9, 3.5517173308154583`*^9}, {
  3.5517283754521756`*^9, 3.5517283755511813`*^9}, {3.552168156124545*^9, 
  3.552168156534569*^9}}],

Cell["Final State S", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.45665974584375*^9, 
  3.456659752203125*^9}, {3.45665980753125*^9, 3.456659829375*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"SFin", "=", 
   RowBox[{"N", "[", 
    FractionBox["1", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665975521875*^9, 3.456659757859375*^9}, {3.4566598095625*^9, 
   3.45665983434375*^9}, {3.536335431601625*^9, 3.536335431726625*^9}, 
   3.53634248861725*^9, {3.5425567824365253`*^9, 3.5425567858841476`*^9}}],

Cell[CellGroupData[{

Cell["Final State J", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659783546875*^9, 
  3.456659787671875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"jFin", "=", 
   RowBox[{"N", "[", 
    FractionBox["3", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.4566597779375*^9, 3.456659780921875*^9}, {3.456661545015625*^9, 
   3.456661545171875*^9}, 3.4569295484212503`*^9, {3.47308863565625*^9, 
   3.473088640640625*^9}, {3.536335439086*^9, 3.536335439226625*^9}, {
   3.542556787927761*^9, 3.542556790751379*^9}, {3.5512819388935547`*^9, 
   3.5512819389873047`*^9}, {3.5517284386257887`*^9, 3.551728438744796*^9}, {
   3.5521681719094477`*^9, 3.5521681723134713`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Final State F Values", "Subsubtitle",
 CellChangeTimes->{{3.45666290321875*^9, 3.45666294725*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"fMinFin", "=", 
   RowBox[{"Abs", "[", 
    RowBox[{"jFin", "-", "i"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fMaxFin", "=", 
   RowBox[{"jFin", "+", "i"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.456662861609375*^9, 3.456662890515625*^9}, {
  3.456662927125*^9, 3.456662957578125*^9}, {3.54906121109375*^9, 
  3.54906122009375*^9}, {3.5490736735625*^9, 3.5490736805625*^9}, {
  3.549247847415039*^9, 3.549247849538086*^9}, {3.5500073128964844`*^9, 
  3.5500073230839844`*^9}, {3.5500073833496094`*^9, 3.5500073842246094`*^9}, {
  3.5512819426591797`*^9, 3.5512819468466797`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"\[Nu]Fin", "=", 
   RowBox[{"SetPrecision", "[", 
    RowBox[{
     RowBox[{"27398.147", " ", "c"}], ",", "20"}], "]"}]}], " ", 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"aFin", "=", 
   RowBox[{"SetPrecision", "[", 
    RowBox[{
     RowBox[{"0.281", " ", 
      SuperscriptBox["10", "6"]}], ",", "20"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"bFin", "=", "0.0"}], ";"}]}], "Input",
 CellChangeTimes->{{3.5512838488154297`*^9, 3.5512838633310547`*^9}, {
   3.551717423621767*^9, 3.551717433431328*^9}, 3.5517275821998043`*^9, 
   3.5517281530804567`*^9, {3.5517283163917975`*^9, 3.551728322063122*^9}, {
   3.55172845647081*^9, 3.5517284682794847`*^9}, {3.552168488942581*^9, 
   3.552168492693796*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Final State Linewidth", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
  3.456659988796875*^9, 3.45665998953125*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Gamma]Fin", "=", 
    FractionBox["1", 
     RowBox[{"285.86", " ", 
      SuperscriptBox["10", 
       RowBox[{"-", "9"}]]}]]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"Seconds", " ", "From", " ", "Heavens"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
  3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
  3.4566600681875*^9, 3.45666007421875*^9}, {3.4566616303125*^9, 
  3.456661631765625*^9}, {3.4569295617493753`*^9, 3.4569295620150003`*^9}, {
  3.45951917796875*^9, 3.45951917865625*^9}, {3.536683236757875*^9, 
  3.536683244726625*^9}, {3.538582219701305*^9, 3.5385822378273416`*^9}, {
  3.5512838421123047`*^9, 3.5512838442841797`*^9}, {3.5517274444239235`*^9, 
  3.5517274454019794`*^9}, {3.551727533652027*^9, 3.5517275372902355`*^9}, {
  3.5521683695107503`*^9, 3.5521683771941895`*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Energies, Matrix Elements and Transition Probablities", "Subtitle",
 CellChangeTimes->{{3.456681533796875*^9, 3.456681554234375*^9}}],

Cell["Formula for Calculating Energies with HFS", "Subsubtitle",
 CellChangeTimes->{{3.45666156509375*^9, 3.45666157378125*^9}, {
  3.456662570765625*^9, 3.4566625805625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Kappa]", "[", 
    RowBox[{"f_", ",", "j_"}], "]"}], "=", 
   RowBox[{
    RowBox[{"f", 
     RowBox[{"(", 
      RowBox[{"f", "+", "1"}], ")"}]}], "-", 
    RowBox[{"j", 
     RowBox[{"(", 
      RowBox[{"j", "+", "1"}], ")"}]}], "-", 
    RowBox[{"i", 
     RowBox[{"(", 
      RowBox[{"i", "+", "1"}], ")"}]}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45666169846875*^9, 3.456661713*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Delta]E", "[", 
    RowBox[{"a_", ",", "b_", ",", "f_", ",", "j_"}], "]"}], ":=", 
   RowBox[{"If", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"j", "\[NotEqual]", " ", 
       FractionBox["1", "2"]}], "&&", 
      RowBox[{"i", "\[NotEqual]", 
       FractionBox["1", "2"]}]}], ",", 
     RowBox[{
      RowBox[{
       FractionBox["1", "2"], "a", " ", 
       RowBox[{"\[Kappa]", "[", 
        RowBox[{"f", ",", "j"}], "]"}]}], "+", 
      RowBox[{"b", " ", 
       FractionBox[
        RowBox[{
         RowBox[{
          FractionBox["3", "4"], 
          SuperscriptBox[
           RowBox[{"\[Kappa]", "[", 
            RowBox[{"f", ",", "j"}], "]"}], "2"]}], "+", 
         RowBox[{
          FractionBox["3", "4"], 
          RowBox[{"\[Kappa]", "[", 
           RowBox[{"f", ",", "j"}], "]"}]}], "-", 
         RowBox[{"i", 
          RowBox[{"(", 
           RowBox[{"i", "+", "1"}], ")"}], "j", 
          RowBox[{"(", 
           RowBox[{"j", "+", "1"}], ")"}]}]}], 
        RowBox[{"2", " ", "i", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", "i"}], "-", "1"}], ")"}], "j", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", "j"}], "-", "1"}], ")"}]}]]}]}], ",", 
     RowBox[{
      FractionBox["1", "2"], "a", " ", 
      RowBox[{"\[Kappa]", "[", 
       RowBox[{"f", ",", "j"}], "]"}]}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.456661820734375*^9, 3.45666183171875*^9}, {
   3.456662393875*^9, 3.456662397109375*^9}, {3.456683422046875*^9, 
   3.456683431546875*^9}, 3.460223934290625*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"\[Nu]GndToInt", "[", 
    RowBox[{"f_", ",", "fPrime_"}], "]"}], "=", 
   RowBox[{"\[Nu]Int", "-", 
    RowBox[{"\[Delta]E", "[", 
     RowBox[{"aGnd", ",", "bGnd", ",", "f", ",", "jGnd"}], "]"}], "+", 
    RowBox[{"\[Delta]E", "[", 
     RowBox[{"aInt", ",", "bInt", ",", "fPrime", ",", "jInt"}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Nu]GndToFin", "[", 
    RowBox[{"f_", ",", "fDPrime_"}], "]"}], "=", 
   RowBox[{"\[Nu]Fin", "-", 
    RowBox[{"\[Delta]E", "[", 
     RowBox[{"aGnd", ",", "bGnd", ",", "f", ",", "jGnd"}], "]"}], "+", 
    RowBox[{"\[Delta]E", "[", 
     RowBox[{"aFin", ",", "bFin", ",", "fDPrime", ",", "jFin"}], "]"}]}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.45686233334375*^9, 3.4568623374375*^9}, 
   3.4572661415462503`*^9, {3.47308885534375*^9, 3.473088860125*^9}, {
   3.4730964853125*^9, 3.473096517875*^9}, {3.473097064765625*^9, 
   3.47309707559375*^9}, {3.536683396664125*^9, 3.53668342911725*^9}, {
   3.5385822740464134`*^9, 3.5385822762455387`*^9}, {3.538582433139513*^9, 
   3.538582435955674*^9}, {3.53858286930646*^9, 3.538583146448312*^9}, {
   3.5399865266491117`*^9, 3.539986527335513*^9}, {3.5425567991442327`*^9, 
   3.5425568137147264`*^9}, {3.5517174795439653`*^9, 
   3.5517175218973875`*^9}, {3.551727944195509*^9, 3.551727944470525*^9}, {
   3.551728158825785*^9, 3.5517281597138357`*^9}, {3.552168184344159*^9, 
   3.5521681859422507`*^9}, 3.552168222350333*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"\[Nu]GndToFin", "[", 
   RowBox[{"2", ",", "3"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.5517281719065332`*^9, 3.551728192390705*^9}}],

Cell[BoxData["8.212026391147076`*^14"], "Output",
 CellChangeTimes->{{3.551728180674035*^9, 3.5517281927097235`*^9}, 
   3.5521685420859375`*^9, 3.552753833303711*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Reduced Matrix Elements in J Basis", "Subsubtitle"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"dRedJGndToInt", "=", "1"}], ";"}], 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     RowBox[{"(", "4.51", ")"}], "1.28", " ", 
     SuperscriptBox["10", "6"]}], " ", 
    RowBox[{"(*", 
     RowBox[{"Hz", "/", 
      RowBox[{"(", 
       RowBox[{"V", "/", "cm"}], ")"}]}], "*)"}], ";"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"dRedJIntToFin", "=", "1"}], ";"}], 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     FractionBox["1", 
      RowBox[{"e", " ", "a0"}]], 
     RowBox[{"\[Sqrt]", 
      RowBox[{"(", 
       RowBox[{
        FractionBox[
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", " ", 
             FractionBox["1", "2"]}], "+", "1"}], " ", ")"}], "3", " ", 
          "hBar", " ", 
          SuperscriptBox["c", "3"]}], 
         RowBox[{"4", " ", 
          SuperscriptBox[
           RowBox[{"(", 
            RowBox[{"2", " ", "\[Pi]", " ", "\[Nu]6P12To8S"}], ")"}], "3"]}]],
         "2.04", " ", 
        SuperscriptBox["10", "6"]}], ")"}]}], "1.28", " ", 
     SuperscriptBox["10", "6"]}], " ", 
    RowBox[{"(*", 
     RowBox[{"Hz", "/", 
      RowBox[{"(", 
       RowBox[{"V", "/", "cm"}], ")"}]}], "*)"}], ";"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.4566056364375*^9, 3.4566056744375*^9}, {
  3.456662686453125*^9, 3.456662724453125*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Reduced Matrix Elements in F Basis", "Subsubtitle",
 CellChangeTimes->{3.4566632656875*^9}],

Cell[BoxData[
 RowBox[{"Do", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Do", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{
       SubscriptBox["dRedGndToInt", 
        RowBox[{"f", ",", "fPrime"}]], "=", 
       RowBox[{"\[Sqrt]", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", " ", "fPrime"}], "+", "1"}], ")"}], 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", "f"}], "+", "1"}], ")"}], 
          SuperscriptBox[
           RowBox[{"SixJSymbol", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{"jInt", ",", "fPrime", ",", "i"}], "}"}], ",", 
             RowBox[{"{", 
              RowBox[{"f", ",", "jGnd", ",", "1"}], "}"}]}], "]"}], "2"], 
          SuperscriptBox["dRedJGndToInt", "2"]}], ")"}]}]}], ";"}], ",", 
     RowBox[{"{", 
      RowBox[{"fPrime", ",", "fMinInt", ",", "fMaxInt"}], "}"}]}], "]"}], ",",
    "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{"f", ",", "fMinGnd", ",", "fMaxGnd"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.456662795671875*^9, 3.456662835546875*^9}, {
  3.456663002453125*^9, 3.456663020171875*^9}, {3.5372063527331247`*^9, 
  3.5372063614518747`*^9}}],

Cell[BoxData[
 RowBox[{"Do", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Do", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      SubscriptBox["dRedIntToFin", 
       RowBox[{"fPrime", ",", "fDPrime"}]], "=", 
      RowBox[{"\[Sqrt]", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", " ", "fDPrime"}], "+", "1"}], ")"}], 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", "fPrime"}], "+", "1"}], ")"}], 
         SuperscriptBox[
          RowBox[{"SixJSymbol", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"jFin", ",", "fDPrime", ",", "i"}], "}"}], ",", 
            RowBox[{"{", 
             RowBox[{"fPrime", ",", "jInt", ",", "1"}], "}"}]}], "]"}], "2"], 
         SuperscriptBox["dRedJIntToFin", "2"]}], ")"}]}]}], ",", 
     RowBox[{"{", 
      RowBox[{"fDPrime", ",", "fMinFin", ",", "fMaxFin"}], "}"}]}], "]"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{"fPrime", ",", "fMinInt", ",", "fMaxInt"}], "}"}]}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.456663051265625*^9, 3.456663093375*^9}, {
   3.456663125046875*^9, 3.456663154296875*^9}, 3.456663192609375*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Matrix Elements in F Basis", "Subsubtitle",
 CellChangeTimes->{{3.456663276328125*^9, 3.456663278140625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Do", "[", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Do", "[", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Do", "[", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"Do", "[", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{
           SubscriptBox["dGndToInt", 
            RowBox[{"f", ",", "fPrime", ",", "m", ",", "mPrime"}]], "=", 
           RowBox[{
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"-", "1"}], ")"}], 
             RowBox[{"fPrime", "-", "mPrime"}]], 
            RowBox[{
             UnderoverscriptBox["\[Sum]", 
              RowBox[{"n", "=", "1"}], "3"], 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{
                "q1", "\[LeftDoubleBracket]", "n", "\[RightDoubleBracket]"}], 
                RowBox[{"ThreeJSymbol", "[", 
                 RowBox[{
                  RowBox[{"{", 
                   RowBox[{"fPrime", ",", 
                    RowBox[{"-", "mPrime"}]}], "}"}], ",", 
                  RowBox[{"{", 
                   RowBox[{"1", ",", 
                    RowBox[{"n", "-", "2"}]}], "}"}], ",", 
                  RowBox[{"{", 
                   RowBox[{"f", ",", "m"}], "}"}]}], "]"}]}], ")"}], 
              SubscriptBox["dRedGndToInt", 
               RowBox[{"f", ",", "fPrime"}]]}]}]}]}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"{", 
           RowBox[{"m", ",", 
            RowBox[{"-", "f"}], ",", "f"}], "}"}]}], "]"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"{", 
         RowBox[{"mPrime", ",", 
          RowBox[{"-", "fPrime"}], ",", "fPrime"}], "}"}]}], "]"}], ",", 
      "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{"fPrime", ",", "fMinInt", ",", "fMaxInt"}], "}"}]}], "]"}], 
    ",", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{"f", ",", "fMinGnd", ",", "fMaxGnd"}], "}"}]}], "]"}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.45666329196875*^9, 3.4566633891875*^9}, {3.456663993*^9,
   3.45666408128125*^9}, {3.45666411909375*^9, 3.456664121515625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Do", "[", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Do", "[", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Do", "[", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"Do", "[", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{
           SubscriptBox["dIntToFin", 
            RowBox[{
            "fPrime", ",", "fDPrime", ",", "mPrime", ",", "mDPrime"}]], "=", 
           RowBox[{
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"-", "1"}], ")"}], 
             RowBox[{"fDPrime", "-", "mDPrime"}]], 
            RowBox[{
             UnderoverscriptBox["\[Sum]", 
              RowBox[{"n", "=", "1"}], "3"], 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{
                "q2", "\[LeftDoubleBracket]", "n", "\[RightDoubleBracket]"}], 
                RowBox[{"ThreeJSymbol", "[", 
                 RowBox[{
                  RowBox[{"{", 
                   RowBox[{"fDPrime", ",", 
                    RowBox[{"-", "mDPrime"}]}], "}"}], ",", 
                  RowBox[{"{", 
                   RowBox[{"1", ",", 
                    RowBox[{"n", "-", "2"}]}], "}"}], ",", 
                  RowBox[{"{", 
                   RowBox[{"fPrime", ",", "mPrime"}], "}"}]}], "]"}]}], ")"}], 
              SubscriptBox["dRedIntToFin", 
               RowBox[{"fPrime", ",", "fDPrime"}]]}]}]}]}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"{", 
           RowBox[{"mPrime", ",", 
            RowBox[{"-", "fPrime"}], ",", "fPrime"}], "}"}]}], "]"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"{", 
         RowBox[{"mDPrime", ",", 
          RowBox[{"-", "fDPrime"}], ",", "fDPrime"}], "}"}]}], "]"}], ",", 
      "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{"fDPrime", ",", "fMinFin", ",", "fMaxFin"}], "}"}]}], "]"}], 
    ",", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{"fPrime", ",", "fMinInt", ",", "fMaxInt"}], "}"}]}], "]"}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.45666329196875*^9, 3.456663489203125*^9}, {
  3.456664105703125*^9, 3.45666415221875*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Doppler Distribution Parameters", "Subtitle",
 FontColor->RGBColor[0, 0.356863, 0]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"amu", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"1.66053886", " ", 
      SuperscriptBox["10", 
       RowBox[{"-", "27"}]]}], ")"}], 
    SuperscriptBox["10", "3"]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"mK", "=", 
   RowBox[{"38.96370668", "amu"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"temp", "=", "343"}], ";"}]}], "Input",
 CellChangeTimes->{{3.456606503203125*^9, 3.4566065051875*^9}, {
   3.456664763890625*^9, 3.456664790515625*^9}, {3.473097352640625*^9, 
   3.473097368078125*^9}, {3.536605457211*^9, 3.53660545749225*^9}, 
   3.5366060920235*^9, {3.536606131101625*^9, 3.536606151414125*^9}, 
   3.53666517880475*^9, {3.540734490795257*^9, 3.5407345006498213`*^9}, {
   3.5407345404901*^9, 3.5407345411041346`*^9}, {3.5517172334698906`*^9, 
   3.5517172625375533`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"\[Alpha]", "=", 
   RowBox[{"\[Sqrt]", 
    RowBox[{"(", 
     FractionBox[
      RowBox[{"2", " ", "kB", " ", "temp"}], "mK"], ")"}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45666479665625*^9, 3.456664796984375*^9}, {
  3.536606154070375*^9, 3.536606154336*^9}, {3.551717264712678*^9, 
  3.5517172690419254`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"dopplerDist", "[", "v_", "]"}], ":=", 
   RowBox[{"Exp", "[", 
    RowBox[{"-", 
     FractionBox[
      SuperscriptBox["v", "2"], 
      SuperscriptBox["\[Alpha]", "2"]]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{3.551717270225993*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Comb Parameters", "Subtitle",
 FontColor->RGBColor[0, 0.356863, 0]],

Cell["Comb Modes", "Subsubtitle",
 CellChangeTimes->{{3.4566813895*^9, 3.45668139175*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"n1NomA", "=", 
   RowBox[{"Round", "[", 
    FractionBox[
     RowBox[{"\[Nu]GndToInt", "[", 
      RowBox[{"fMinGnd", ",", "fMinInt"}], "]"}], "fRepAvg"], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"n1NomB", "=", 
   RowBox[{"Round", "[", 
    FractionBox[
     RowBox[{"\[Nu]GndToInt", "[", 
      RowBox[{"fMaxGnd", ",", "fMaxInt"}], "]"}], "fRepAvg"], "]"}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.456665147765625*^9, 3.456665238296875*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"n2NomA", "=", 
   RowBox[{"Round", "[", 
    RowBox[{
     FractionBox["1", "fRepAvg"], 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"\[Nu]GndToFin", "[", 
        RowBox[{"fMinGnd", ",", "fMinFin"}], "]"}], "-", 
       RowBox[{"\[Nu]GndToInt", "[", 
        RowBox[{"fMinGnd", ",", "fMinInt"}], "]"}]}], ")"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"n2NomB", "=", 
   RowBox[{"Round", "[", 
    RowBox[{
     FractionBox["1", "fRepAvg"], 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"\[Nu]GndToFin", "[", 
        RowBox[{"fMaxGnd", ",", "fMaxFin"}], "]"}], "-", 
       RowBox[{"\[Nu]GndToInt", "[", 
        RowBox[{"fMaxGnd", ",", "fMaxInt"}], "]"}]}], ")"}]}], "]"}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.456665147765625*^9, 3.4566653128125*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"n1List", "=", 
  RowBox[{"Union", "[", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"n1NomA", "+", "\[Delta]n"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Delta]n", ",", 
         RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"n1NomB", "+", "\[Delta]n"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Delta]n", ",", 
         RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "]"}]}], "]"}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.456665383359375*^9, 3.45666546196875*^9}, {
  3.456665618125*^9, 3.456665620609375*^9}, {3.5372066701081247`*^9, 
  3.5372066730768747`*^9}, {3.5377187815151253`*^9, 3.5377187956088753`*^9}, {
  3.5377946040932503`*^9, 3.5377946077338753`*^9}, {3.5377949385776253`*^9, 
  3.5377949450151253`*^9}, {3.537893207650926*^9, 3.537893217432143*^9}, {
  3.5385234618671618`*^9, 3.53852346744048*^9}, {3.538583211301021*^9, 
  3.5385832213905983`*^9}, {3.5407344378562293`*^9, 3.5407344440115814`*^9}, {
  3.5409093630926204`*^9, 3.540909373629223*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "1563369", ",", "1563370", ",", "1563371", ",", "1563372", ",", "1563373", 
   ",", "1563374", ",", "1563375", ",", "1565216", ",", "1565217", ",", 
   "1565218", ",", "1565219", ",", "1565220", ",", "1565221", ",", 
   "1565222"}], "}"}]], "Output",
 CellChangeTimes->{3.5521685429609375`*^9, 3.552753835763672*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"n2List", "=", 
  RowBox[{"Union", "[", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"n2NomA", "+", "\[Delta]n"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Delta]n", ",", 
         RowBox[{"-", "4"}], ",", "4"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"n2NomB", "+", "\[Delta]n"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Delta]n", ",", 
         RowBox[{"-", "4"}], ",", "4"}], "}"}]}], "]"}]}], "]"}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.456665383359375*^9, 3.45666546196875*^9}, {
  3.456665585125*^9, 3.45666558903125*^9}, {3.456665624328125*^9, 
  3.456665626734375*^9}, {3.5372066753268747`*^9, 3.5372066804518747`*^9}, {
  3.5377187866401253`*^9, 3.5377187978745003`*^9}, {3.5377946088901253`*^9, 
  3.5377946116245003`*^9}, {3.5377949469370003`*^9, 3.5377949552338753`*^9}, {
  3.5378932198813477`*^9, 3.537893225029357*^9}, {3.538523470201638*^9, 
  3.5385234737538414`*^9}, {3.538583224257762*^9, 3.5385832272179317`*^9}, {
  3.540734445712679*^9, 3.540734450881974*^9}, {3.540909375727343*^9, 
  3.5409093833457785`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "1721434", ",", "1721435", ",", "1721436", ",", "1721437", ",", "1721438", 
   ",", "1721439", ",", "1721440", ",", "1721441", ",", "1721442"}], 
  "}"}]], "Output",
 CellChangeTimes->{3.5521685436171875`*^9, 3.552753836354492*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Comb Equation", "Subsubtitle",
 CellChangeTimes->{{3.456681410328125*^9, 3.456681413015625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Omega]L", "[", 
    RowBox[{"n_", ",", "fR_", ",", "f0_"}], "]"}], ":=", 
   RowBox[{"2", " ", "\[Pi]", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"n", " ", "fR"}], "+", "f0"}], ")"}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.456678515171875*^9, 3.45667856759375*^9}, {
   3.4566795630625*^9, 3.45667956609375*^9}, 3.5413431977628536`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Integrate 2-photon Lineshape over Doppler Distribution", "Subtitle",
 FontColor->RGBColor[0, 0.356863, 0]],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"transProb", "[", 
     RowBox[{"\[Omega]1_", ",", "\[Omega]2_", ",", "v_"}], "]"}], "=", 
    RowBox[{
     UnderoverscriptBox["\[Sum]", 
      RowBox[{"fDPrime", "=", "fMinFin"}], "fMaxFin"], 
     RowBox[{
      UnderoverscriptBox["\[Sum]", 
       RowBox[{"f", "=", "fMinGnd"}], "fMaxGnd"], 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          FractionBox[
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"\[Gamma]Fin", "+", "\[Gamma]T"}], ")"}], "2"], "4"], 
          "/", 
          RowBox[{"(", 
           RowBox[{
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{
               RowBox[{"2", " ", "\[Pi]", " ", 
                RowBox[{"\[Nu]GndToFin", "[", 
                 RowBox[{"f", ",", "fDPrime"}], "]"}]}], "-", 
               RowBox[{"\[Omega]1", 
                RowBox[{"(", 
                 RowBox[{"1", "+", 
                  FractionBox["v", "c"]}], ")"}]}], "-", 
               RowBox[{"\[Omega]2", 
                RowBox[{"(", 
                 RowBox[{"1", "-", 
                  FractionBox["v", "c"]}], ")"}]}]}], ")"}], "2"], "+", 
            FractionBox[
             SuperscriptBox[
              RowBox[{"(", 
               RowBox[{"\[Gamma]Fin", "+", "\[Gamma]T"}], ")"}], "2"], "4"]}],
            ")"}]}], ")"}], 
        FractionBox["1", 
         RowBox[{
          RowBox[{"2", "f"}], "+", "1"}]], 
        RowBox[{
         UnderoverscriptBox["\[Sum]", 
          RowBox[{"m", "=", 
           RowBox[{"-", "f"}]}], "f"], 
         RowBox[{
          UnderoverscriptBox["\[Sum]", 
           RowBox[{"mDPrime", "=", 
            RowBox[{"-", "fDPrime"}]}], "fDPrime"], 
          SuperscriptBox[
           RowBox[{"Abs", "[", 
            RowBox[{
             UnderoverscriptBox["\[Sum]", 
              RowBox[{"fPrime", "=", "fMinInt"}], "fMaxInt"], 
             RowBox[{
              UnderoverscriptBox["\[Sum]", 
               RowBox[{"mPrime", "=", 
                RowBox[{"-", "fPrime"}]}], "fPrime"], "\[IndentingNewLine]", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{
                 FractionBox[
                  RowBox[{"\[Gamma]Int", "+", "\[Gamma]T"}], "2"], "/", 
                 RowBox[{"(", 
                  RowBox[{
                   FractionBox[
                    RowBox[{"\[ImaginaryI]", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Gamma]Int", "+", "\[Gamma]T"}], ")"}]}], "2"], 
                   "+", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"\[Omega]1", 
                    RowBox[{"(", 
                    RowBox[{"1", "+", 
                    FractionBox["v", "c"]}], ")"}]}], "-", " ", 
                    RowBox[{"2", " ", "\[Pi]", " ", 
                    RowBox[{"\[Nu]GndToInt", "[", 
                    RowBox[{"f", ",", "fPrime"}], "]"}]}]}], ")"}]}], ")"}]}],
                 ")"}], 
               SubscriptBox["dGndToInt", 
                RowBox[{"f", ",", "fPrime", ",", "m", ",", "mPrime"}]], 
               SubscriptBox["dIntToFin", 
                RowBox[{
                "fPrime", ",", "fDPrime", ",", "mPrime", ",", 
                 "mDPrime"}]]}]}]}], "]"}], "2"]}]}]}], ")"}]}]}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"fnToIntegrate", "[", 
     RowBox[{"fR_", ",", "fOffset_", ",", "v_"}], "]"}], "=", 
    RowBox[{
     UnderoverscriptBox["\[Sum]", 
      RowBox[{"n1", "=", "1"}], 
      RowBox[{"Length", "[", "n1List", "]"}]], 
     RowBox[{
      UnderoverscriptBox["\[Sum]", 
       RowBox[{"n2", "=", "1"}], 
       RowBox[{"Length", "[", "n2List", "]"}]], 
      RowBox[{
       RowBox[{"dopplerDist", "[", "v", "]"}], 
       RowBox[{"transProb", "[", 
        RowBox[{
         RowBox[{"\[Omega]L", "[", 
          RowBox[{
           RowBox[{
           "n1List", "\[LeftDoubleBracket]", "n1", "\[RightDoubleBracket]"}], 
           ",", "fR", ",", "fOffset"}], "]"}], ",", 
         RowBox[{"\[Omega]L", "[", 
          RowBox[{
           RowBox[{
           "n2List", "\[LeftDoubleBracket]", "n2", "\[RightDoubleBracket]"}], 
           ",", "fR", ",", "fOffset"}], "]"}], ",", "v"}], "]"}]}]}]}]}], 
   ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"vIntSpec", "[", 
    RowBox[{
     RowBox[{"fRep_", "?", "NumericQ"}], ",", 
     RowBox[{"fO_", "?", "NumericQ"}]}], "]"}], ":=", 
   RowBox[{"NIntegrate", "[", 
    RowBox[{
     RowBox[{"fnToIntegrate", "[", 
      RowBox[{"fRep", ",", "fO", ",", "v"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"v", ",", 
       RowBox[{"-", "60000"}], ",", "60000"}], "}"}], ",", 
     RowBox[{"MinRecursion", "\[Rule]", "3"}], ",", 
     RowBox[{"MaxRecursion", "\[Rule]", "10"}], ",", 
     RowBox[{"AccuracyGoal", "\[Rule]", "4"}], ",", 
     RowBox[{"PrecisionGoal", "\[Rule]", "2"}]}], "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.538523548469115*^9, 3.5385235524603434`*^9}, {
   3.5385832601238136`*^9, 3.538583261230877*^9}, {3.540909398014618*^9, 
   3.540909411994417*^9}, 3.541343202952861*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Calculation", "Subtitle",
 CellChangeTimes->{{3.456681446140625*^9, 3.45668144753125*^9}, {
  3.5379007941050625`*^9, 3.537900798020669*^9}},
 FontColor->RGBColor[0, 0.356863, 0]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"rFrep", "=", "0.0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"tt", "=", 
    RowBox[{"AbsoluteTime", "[", "]"}]}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{
   "distribute", " ", "NIntegrate", " ", "among", " ", "CPU", " ", "cores"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"specList", "=", 
   RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"SetSharedVariable", "[", "specList", "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ParallelDo", "[", 
   RowBox[{
    RowBox[{"AppendTo", "[", 
     RowBox[{"specList", ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{
          "fRepList", "\[LeftDoubleBracket]", "kk", "\[RightDoubleBracket]"}],
           "-", "rFrep"}], ")"}], ",", 
        RowBox[{"Chop", "[", 
         RowBox[{"vIntSpec", "[", 
          RowBox[{
           RowBox[{
           "fRepList", "\[LeftDoubleBracket]", "kk", 
            "\[RightDoubleBracket]"}], ",", 
           RowBox[{
           "fOffsetList", "\[LeftDoubleBracket]", "kk", 
            "\[RightDoubleBracket]"}]}], "]"}], "]"}]}], "}"}]}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"kk", ",", 
      RowBox[{"Length", "[", "fRepList", "]"}]}], "}"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ttt", "=", 
   RowBox[{
    RowBox[{"AbsoluteTime", "[", "]"}], "-", "tt"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Print", "[", 
  RowBox[{"\"\<Execution time = \>\"", ",", "ttt"}], "]"}]}], "Input",
 CellChangeTimes->{{3.5407346692514644`*^9, 3.5407346703575277`*^9}, {
  3.5413432101228704`*^9, 3.5413432306928997`*^9}, {3.5425568547897897`*^9, 
  3.542556855210992*^9}, {3.5492478888847656`*^9, 3.549247903041992*^9}, {
  3.552168231367849*^9, 3.5521682316408644`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Execution time = \"\>", "\[InvisibleSpace]", 
   "1689.265625`10.679242938153749"}],
  SequenceForm["Execution time = ", 1689.265625`10.679242938153749],
  Editable->False]], "Print",
 CellChangeTimes->{3.552170227366211*^9, 3.5527555276621094`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{
   RowBox[{"Sort", "[", "specList", "]"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
   RowBox[{"Joined", "\[Rule]", "True"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.5407062102656565`*^9, 3.5407062108584414`*^9}, 
   3.5407346912477226`*^9, {3.5407454061664057`*^9, 3.54074540677481*^9}, 
   3.540917370755953*^9, {3.540955271524134*^9, 3.5409552734745717`*^9}, {
   3.5409573772453723`*^9, 3.5409573788055224`*^9}, 3.5413432386329107`*^9}],

Cell[BoxData[
 GraphicsBox[
  {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwl2HdcT/sDx/HTQEoppTLSIg3lkpJKJQkJRZQUQiipVEKKo0GpKE0z1yYZ
2SujrIy6RVFGSPPHDZE07u++39c/Hs/H69v3++mczzjftBcHzfKVFATBTkIQ
8D/++T7NmyA3Qz3A7Xx/O0EQ6YChHrs/7Fb914IEPKTv8aDLZ2GRjl5TsClm
CV4vSMLtzf0NfIP6odN29j2EjO3y6FLwM0+39ctNeqHTTVXbs1r3/Lb9t0vD
aiZHMn/Gdv9rka7Url8Z3N2O3oPvN3mtepJ0GzqdfE+nqWjmD/Se8IowsweP
YlvR6cjiWkV7VfZecGzvnB3rNv5Epx1fjHbot/sXugxcYlwctDkZ4xHpWfOV
TMwtO9B7w2OfL1FoHsBOS1p3XNr/J8YjyMLjPFZWDyvh59MRuy5O2ZL2DV0O
flilpysZ9BWdzih/XxOzFhb6wObZgyrsjuH1Im3cpfKsSIbjl4e7+63bfngl
xivSr3vtTGuv6ERXgH/USnT7nBdwfemLYXt7lsVI4vr3hZePtjozapkUOl0z
p+f29BxYUIQXTHBLsG7C60X6kIZf9Y8CvJ+gxPsVnhflvgDXQ6QlalpTLVZy
fP3gvVuW9rt4qwWdluiTWePyrAldGR5+e/r02m+f0Ol5Rckvy4I/oKvAp/q7
Tvq45R06/W3D01UabdXo/eFhw+5dVpB4hU4vMtXTyMivQFfl/bLUEbUKn6PT
CbnnZ4QllKOrwTKHWnbIGZah02uiU29MNPkLXR3OC5t7zW9DKTp9e77+waFq
sDAAlm+7lPoguwSdLg01OjYtFxYGwp72tiolQ/nz9MYeWWYf3vPnB8HW1de6
/niDzxPpgKihD8P+wXiEwbDH5+kuK5swXpFOuJK84Ukhx68B99f5oVYdzvHT
1skXhpwM4vsPgevOWV+dpPIUnQ4NDpPRP/wIXRPWjLp+ytXmATrtrz7cVyfq
HroWx+Pw4lvpwUJ0eun1D6u2fr2Lrg2bPbz69+lpsEgbPni5b+ChO+g6XK8b
2o+XKMMifepZmc3iPbfRdeFftZa3LSRgkd7UZuCme+0W+lC+/6jKeUcHwCId
kmlxZerMAvRhsIllt2SmFCzS02LaU/ck3UTXg/8XkBF3QxsW6S8Pa+Ibm66j
D+f1qXvkUHEUFukjflMnFXaw68NbI6TTD5xmp5tHmGifCWM3gG8m1krMNWOn
r8s/qMvox24Ij55l3BHWcA2dfmQcqO9bCAtG8MebPesn57DTG9cqjev0Zh8B
25QnZ831ZafV76XoNBuxG8Mvv0sNe67FTs+2HHXQaiK7CayjbdquZ8dOL/Q7
Fuk9nn0k50NMzu7CCey0YsWnbenO7H9w/j3elnmM4xHpp81VC93WsI/i/pmg
3FmYyk5vXhq8I2Y9+2i43sc4OK6CnR53MTn45BReH1P4VfbnlPUnef3o3x+O
L7S/wT6G+5v6OrUDe9npaPedOUoiuxk/f11rT2+9G+h0hq2eU5IHLJjDmW9D
+jrks9MnD9Qrak7n/Bn73/7gOyV6P+cH7W1wYfwNbXYL+Iyv9QubHZxfdNUu
p9iaL+zj4BlKYz8P28/5SVsvMa5sHM35bcn5fMr61YA8zm86+ODjKt0ydivu
Dz7eAyVnc33Qmx1Hdh8exvVlDdt+d3olXcL1Rc8YMtBpUBLWozAeVnrdcqev
RxE63bhVN6fXaq5fG1j/vtfrU6fvo9Npy+cvlrN+iG4Lqy5OuvtmXDE6/bv8
WrDL/Sfodrx+qklpfhbcv+i5X6PqM7K4P03geVF/Of+sLvZbkZaJ8Kh7V1mF
bg/fCbpRl2D4Hp2e4yU1o9AH+70wEXaRPTrDpbEBnT6S+6KjurUR3YH31/WK
iZIyOz1tslSVdmQt+iSupyE+wy/9wfenc6d6b1/+9A26I5zu7/DzlBPGI9Ix
dhPNrlq9RJ8MO6W0JDaHcfx01UPH4+vtX6BPgWuloy8HSsMi3ejwvnG8B/tU
WPplmdXBJna6T0eOq24Kzx8n2Erd/5fs35XotPatj88NN+K8EqbBz8PDrQ0e
8PyiBzc6vlDSeovuzPWy9PXbZRdr0P/zmHPObl94Pk7n+fvJxvLvLzw/6eRN
N47vquD1mwEPaast2lf3BZ0WZHde0tf7jj4T9nmcLbtoLp8faPWREfFpaf+g
u3C+NDk3u6yWxvlPV79cUmkW0Rvnvyuvt2nXvCJ3PM+J9Gn/vy+MiFNAnwVv
M9rtt/FsH3T6Zr9j+xI/yKDPhj/fsspqzOL70xYVjpVKrXy+cON8sHL8UmCO
5xmRLnD7+L+WcXz+mwOPCnuratvF5z9a9Xe7T14rn7/mcj0+drh9Nou/H+0e
WzIpYQV/3p3zu0v7ULosn9/oGukqD0M/Pn96wEVFLSqH5dlph3ENWedS+Xw4
j/vtIcP9vp/4/EhvyZzt+OEen388+flunckbwmGRHnS4c3mSEl4vzIcnFldJ
K+7h+OkTzVHfDTz5+V7cH/O9Ol5b8vmXfp0yOyKrQgLXxxs2KM4f1FXZA9eP
znjcR80ySQ59Ae/nnwn78v5SQqeP6/6KNV2ohr6Q83td9Nhekhro9AOrhmoN
Ly30RfDAkfMXiF9hkY7JC3e2z9NE9+H9GhwfqDZzEDqtFuLufW0Xvx8s5v7n
1xVn5MrnfzrnyLXri8L6oi/h+Ry/M7MljvODtvi8/GHHUY5/KSyRcL7zoy07
3cN5/KH8Gn5/8IWz+vp8NI3D+4l0/XXNVd7rFdGXcT2uXauiWgyL9L71bs5r
ktiX8/Vfnuy9a435KtK2fde6bDHj56+AzS3VC5zf8/sJ3eo/tyUxBddb8IOz
V9i8DXzH+Uu7dTyaX/sNFvx5/WzntUjp9kSnO5JrSn8v4vxfCXdaX4r/KPD3
o9tcR0aU3+T4AuA4z/r2RwtV0Ol9vkZHwqUGoq+CHSNCWvzP8f7QJ+S1DI3D
hqIHwvGBt7fvnz0cnV5zNy5foxsWgrieFmeX2LvqodPT3VcNNE/URQ+G15lf
um+kz/eng39bhcx+ys9fzevv+K37zCTMJ5HW2XrhfuhRfv8L4X7wtMho/il+
f6SdNeruBnRzfoTCJl8dH+frDECnj5vu+h6/CPNJCON5fOZxYoXA+UlH27dI
dTfDwhrurwn67hGf2GnJ3Y8CLisMRg/nekhcOFQmVB2dDvHuTpwWj+sprOV5
kpT9RsGT64O22aWVkp7K/WsdvLT0mYGSK+8PnXmiNuVLLufHevjTs4PtNgXs
dNHSUXZ7NTm/I+DQuy2ZUrf4/rSff4z/bX1enw3wFZMNO+O1cD1Feozzm6pi
Ta6/SFg5VTbmUTruh0gfNymbrKSM+yVEwV6Dzn2STOP9o2etTXh5ZcQw9I2w
fGJy28YFOuj0jil9LO7YD0HfBD9x8J+iLcfPp41y7C3C7bk/iFy/sre2nrfn
93s6bGRxgJ0K+2Y4weuOYrAT7x9ddlmz9uYWXv9ouNwj95VTK/YLkXaYY1E3
1V8bPQY++nbEtpHGHB8tv8VZMXoJLMTCUX76l7v34vUi/dem+6e+SXD/iYO/
GoSOnRjEz6d72By7KxvI+7uF98uydeXdV1z/9EUrhSeHFbl/bOV+c+xAbbG1
LDrt2VVq0/GD6zMejmxY2BXaDIt0Z45C4IMovF5IgBeYhMj8us/zj+5dMLUw
OgT3W9jG/SvSwFFNhvOfPt3hvaT3W8xHIRE+ZNS8p3Eb9086NrKlvKuU8z8J
toyz89P/zN+Pbli1oGn8S66fZNjMI7RtyTJldHqr16j66SO4f2yHPTbG5ZXK
cX+kv/9u23U/jfN7B7wiWEWjeiKvD31Hp+xIgzP2ayEF/qiQn9FqzPHTg7eN
GD3Gg+NPhSutS8c+kef46d4Du1XjT/D+7+T++87mZ4Q71yc9evPwBt3z7Gk8
39883d/7CX8/WnmzV3kfT66PdI4v1cXqagnHR68xDIn2LuT1z+D8cL+uNaWc
5x89smGAeWCsJHomLH0g/IZeAP/+QWf7mLuop+E8FbLgW+oqz8/87ML5SpvW
DV6d8Q97Nsc/XCEyNZU/T+fsGF0a5sT33wUvvmLa3zObn08/17LvCpiE80LY
zfWVdsxonBHPD7pX/qefPar595U93C+UezyRm8e/v9CJGRKFh3/z+WcvHDg8
t+emaTjvRTq3f8zJ5Z38+9M++HWZ+nulORw/rbwyuj2ggM9v+7nepgxNS3/O
n6czXFXTI3/xfMrhfDHTmLk0nfObPhiV227xXs7u/59g/tU=
    "]]},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  PlotRangeClipping->True]], "Output",
 CellChangeTimes->{3.552170227397461*^9, 3.5527555277685547`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Export", "[", 
  RowBox[{"oFile", ",", 
   RowBox[{"Sort", "[", "specList", "]"}], ",", "\"\<Table\>\""}], 
  "]"}]], "Input"],

Cell[BoxData["\<\"C:\\\\Users\\\\work\\\\Data\\\\output\\\\K39-D2-4D32.txt\"\>\
"], "Output",
 CellChangeTimes->{3.552170227694336*^9, 3.552175407616211*^9, 
  3.552755528203125*^9}]
}, Open  ]]
}, Open  ]]
},
WindowSize->{1272, 684},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
FrontEndVersion->"8.0 for Microsoft Windows (32-bit) (November 7, 2010)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[545, 20, 486, 13, 72, "Input"],
Cell[CellGroupData[{
Cell[1056, 37, 407, 9, 85, "Input"],
Cell[1466, 48, 350, 9, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[1853, 62, 132, 2, 28, "Subsubtitle"],
Cell[1988, 66, 1966, 32, 31, "Input"],
Cell[3957, 100, 657, 18, 29, "Text"],
Cell[4617, 120, 4468, 92, 332, "Input"],
Cell[9088, 214, 660, 19, 29, "Text"],
Cell[9751, 235, 2672, 75, 314, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[12460, 315, 76, 1, 49, "Subtitle"],
Cell[12539, 318, 37, 0, 28, "Subsubtitle"],
Cell[12579, 320, 191, 6, 31, "Input"],
Cell[CellGroupData[{
Cell[12795, 330, 40, 0, 28, "Subsubtitle"],
Cell[12838, 332, 223, 8, 31, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[13098, 345, 107, 1, 28, "Subsubtitle"],
Cell[13208, 348, 330, 10, 31, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[13587, 364, 75, 1, 49, "Subtitle"],
Cell[13665, 367, 102, 1, 49, "Subtitle"],
Cell[13770, 370, 149, 2, 28, "Subsubtitle"],
Cell[13922, 374, 357, 8, 47, "Input"],
Cell[CellGroupData[{
Cell[14304, 386, 105, 1, 28, "Subsubtitle"],
Cell[14412, 389, 406, 10, 47, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[14855, 404, 101, 1, 28, "Subsubtitle"],
Cell[14959, 407, 219, 6, 31, "Input"],
Cell[15181, 415, 167, 5, 31, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[15397, 426, 106, 1, 49, "Subtitle"],
Cell[15506, 429, 222, 3, 28, "Subsubtitle"],
Cell[15731, 434, 309, 7, 47, "Input"],
Cell[CellGroupData[{
Cell[16065, 445, 199, 3, 28, "Subsubtitle"],
Cell[16267, 450, 285, 7, 47, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[16589, 462, 109, 1, 28, "Subsubtitle"],
Cell[16701, 465, 630, 13, 52, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[17368, 483, 389, 5, 28, "Subsubtitle"],
Cell[17760, 490, 1129, 20, 31, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[18926, 515, 413, 6, 28, "Subsubtitle"],
Cell[19342, 523, 589, 10, 31, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[19980, 539, 356, 10, 52, "Subtitle"],
Cell[20339, 551, 249, 3, 28, "Subsubtitle"],
Cell[20591, 556, 328, 7, 47, "Input"],
Cell[CellGroupData[{
Cell[20944, 567, 199, 3, 28, "Subsubtitle"],
Cell[21146, 572, 463, 9, 47, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[21646, 586, 115, 1, 28, "Subsubtitle"],
Cell[21764, 589, 347, 9, 52, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[22148, 603, 252, 3, 28, "Subsubtitle"],
Cell[22403, 608, 768, 18, 47, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[23208, 631, 315, 4, 28, "Subsubtitle"],
Cell[23526, 637, 1307, 24, 54, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[24870, 666, 350, 5, 28, "Subsubtitle"],
Cell[25223, 673, 1150, 21, 33, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[26410, 699, 371, 5, 28, "Subsubtitle"],
Cell[26784, 706, 1099, 20, 33, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[27932, 732, 493, 12, 52, "Subtitle"],
Cell[28428, 746, 240, 3, 28, "Subsubtitle"],
Cell[28671, 751, 403, 8, 47, "Input"],
Cell[CellGroupData[{
Cell[29099, 763, 199, 3, 28, "Subsubtitle"],
Cell[29301, 768, 606, 11, 47, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[29944, 784, 104, 1, 28, "Subsubtitle"],
Cell[30051, 787, 641, 13, 52, "Input"],
Cell[30695, 802, 777, 20, 75, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[31509, 827, 293, 4, 28, "Subsubtitle"],
Cell[31805, 833, 973, 19, 47, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[32827, 858, 139, 1, 49, "Subtitle"],
Cell[32969, 861, 175, 2, 28, "Subsubtitle"],
Cell[33147, 865, 442, 15, 31, "Input"],
Cell[33592, 882, 1602, 47, 56, "Input"],
Cell[35197, 931, 1508, 31, 52, "Input"],
Cell[CellGroupData[{
Cell[36730, 966, 194, 4, 52, "Input"],
Cell[36927, 972, 167, 2, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[37131, 979, 57, 0, 28, "Subsubtitle"],
Cell[37191, 981, 1384, 45, 76, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[38612, 1031, 97, 1, 28, "Subsubtitle"],
Cell[38712, 1034, 1254, 33, 99, "Input"],
Cell[39969, 1069, 1184, 32, 99, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[41190, 1106, 115, 1, 28, "Subsubtitle"],
Cell[41308, 1109, 2134, 54, 275, "Input"],
Cell[43445, 1165, 2166, 55, 275, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[45660, 1226, 89, 1, 49, "Subtitle"],
Cell[45752, 1229, 856, 20, 74, "Input"],
Cell[46611, 1251, 352, 9, 53, "Input"],
Cell[46966, 1262, 288, 9, 50, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[47291, 1276, 73, 1, 49, "Subtitle"],
Cell[47367, 1279, 90, 1, 28, "Subsubtitle"],
Cell[47460, 1282, 507, 15, 87, "Input"],
Cell[47970, 1299, 833, 25, 87, "Input"],
Cell[CellGroupData[{
Cell[48828, 1328, 1121, 25, 31, "Input"],
Cell[49952, 1355, 358, 7, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[50347, 1367, 1168, 26, 31, "Input"],
Cell[51518, 1395, 273, 6, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[51828, 1406, 102, 1, 28, "Subsubtitle"],
Cell[51933, 1409, 394, 10, 31, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[52376, 1425, 112, 1, 49, "Subtitle"],
Cell[52491, 1428, 5254, 138, 257, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[57782, 1571, 185, 3, 49, "Subtitle"],
Cell[CellGroupData[{
Cell[57992, 1578, 1856, 52, 152, "Input"],
Cell[59851, 1632, 298, 6, 23, "Print"]
}, Open  ]],
Cell[CellGroupData[{
Cell[60186, 1643, 508, 9, 31, "Input"],
Cell[60697, 1654, 4587, 78, 234, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[65321, 1737, 150, 4, 31, "Input"],
Cell[65474, 1743, 182, 3, 30, "Output"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
